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Abstract 

The notion of the integral over anticommuting Grassmann variables (non- 
quantum fermionic fields) seems to be the most powerful tool in order to 
extract the exact solutions for the 2D Ising models on simple and more com- 
plicated lattices, which is the subject of a discussion in this report. 

1. Introduction. The fermionic structures in the two-dimensional (2D) Ising 
model [1] were first recognized within the transfer-matrix and combinatorial ap- 
proaches [2] . It was realized later on that the notion of the integral over anticom- 
muting Grassmann variables due to Bcrczin [3] is a powerful tool to study the 2D 
Ising model [3-8]. A simple method of fermionic analysis for the 2D Ising models 
disposed on simple and more complicated lattices has been developed in [7-8]. The 
approach is based on the integration over the anticommuting Grassmann variables 
(nonquantum fermionic fields) and the mirror-ordering factorization principle for 
the 2DIM density matrix. The method is straightforward, the traditional transfer- 
matrix or combinatorial considerations are not used. Schematically, we have: 

Q = Sp (5 (cr) — »■ Sp Q (cr I a) ^ Sp (5 (a) = Q ■ (^ \ 

(a) {a\a) (a) 

We start here with the original Ising spin partition function, Q, and introduce, in a 
special way, a set of new purely anticommuting Grassmann variables (a) thus pass- 
ing to a mixed (c | a) representation. Eliminating spin variables (cr) in this mixed 
representation, we obtain a purely fermionic expression for the same partition func- 
tion Q. The final expression for Q appears as a fermionic Gaussian integral. In 
essence, this means the exact solution of the problem. The partition function is 
expressible here as a fermionic Gaussian integral even for the most general inho- 
mogeneous distribution of coupling parameters [7]. In particular, for the standard 
rectangular lattice this gives a few line derivation of the Onsager result [7]. The 
2D Ising lattices with more complicated local structures have been analyzed by the 
factorization method as well [8,9]. The Gaussian fermionic representations have re- 
cently been constructed also for the inhomogeneous 2D dimer models [10]. In this 
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scheme, we introduce Grassmann variables in order to decouple the bond Boltzmann 
weights into separable factors called the Grassmann factors (GFs). We then com- 
bine GFs with the same spin variables into separable groups and sum over spins in 
each group independently thus passing to a purely fermionic representation for Q. 
The key point of all the construction is the mirror-ordering procedure for the global 
products of GFs which enables one to perform the actual summing over spin degrees 
of freedom [7,8] . In what follows we comment shortly on some aspects of the method. 

2. Grassmann vairiables. Grassmann variables are the purely anticommuting 
fermionic numbers. Given a set of Grassmann variables ai, a2, ap^, we have aittj + 
ajtti — 0, Oj — 0. Berezin's rules of integration for one variable are [3]: 

/ daj ■ aj = 1 , J daj -1 = 0. (2) 

In the multidimensional integral, the differentials dai, da2, da^ are again anticom- 
muting with each other and with the variables. The basic relations of the grassman- 
nian analysis concern the Gaussian fermionic integrals [3] . The Gaussian integral of 
the first kind is related to the determinant: 

N f N N \ 

J n da*daj exp J2 J2 cuAija* — det A , (3) 
j=i \i=i j=i J 

where {aj,a*} is a set of completely anticommuting Grassmann variables, the ma- 
trix in the exponential is arbitrary. The fermionic exponential here is assumed in 
the sense of its series expansion, the series terminates at some stage due to the prop- 
erty ttj — 0. The origin of the determinant in (3) is due to the known interrelation 
between the fermionic algebra and determinant combinatorics. By convention, the 
variables and a* can be considered as complex conjugated fermionic fields, oth- 
erwise these are simply independent variables. The Gaussian integral of the second 
kind, for real fermionic fields, is related to the Pfafhan: 

f N N \ 

J dttN ■■■ da2dai exp , aiAijUj — PfaffA, Aij — —Aji, (4) 

V i=ij=i J 

where the matrix is assumed to be skew-symmetric. The Pfaffian is some combina- 
torial polynomial in elements A^^j known in mathematics for a long time. Pfaffian 
combinatorics is in fact identical with that of the fermionic version of Wick's theo- 
rem. For any skew-symmetric matrix {A^j = —Aji) we have: det A = (PfaffA)^ . 

3. The 2D Ising model on rectangular lattice. To illustrate (1), let us 
consider the 2D Ising model on a rectangular lattice net with the inhomogeneous 
distribution of coupling parameters [7]. The hamiltonian is given as follows: 

—(3H(a)— J2 J2 [ ^m+lnC"mnCm+ln + ^mn+l'^mn'^mn+l ] , P—l/kT, (5) 
m=l n=l 

where = J^^/kT, are the magnetic exchange energies, P — 1/kT \s the 
inverse temperature. The Ising spins (7^„ = ±1 are disposed at lattice sites mn, with 
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■m,n — 1, ...,L, N — ^ oo. For finite N, we assume free-boundary conditions: 
CM+in — CmAT+i = 0. Noting that for typical bond weight e^"'"'' — coshft + sinhftcrcr', 
which readily follows from (cro"')^ = 1, we come to the reduced partition function: 

Q = Sp { n ( 1 + ^L+ln'^mnC"m+ln) (1 + t^J^+l'^mnO'mn+l) } , (Q) 
{a) ^ ^ 

where t^^l = tanh 6^^, and Sp(^a) = 2~^S(o-) stands for a properly normalized spin 
averaging [7]. Following (1), we introduce a set of Grassmann variables {amn, 
bmn, ^mn}' ^ P^^^ P®^ bond, and factorize the local bond weights as follows: 

1 "I" t\ran+l^mn(^mn+l — J dh^^dhmn^^"^"^"^'^ "I" ^mn^mn) "I" ^mn^mn+l) ■ 

These identities can be checked simply by using the rules (2), and noting that e''" = 
1 + aa*, since (aa*)^ = 0. Neglecting the sign of the Gaussian fermionic averaging, 
we see that the bond weights are presented now as ^rnn^m+in ' ^mn^mn+n where 
the Grassmann factors , , B^^ , B^^^^ are to be identified ^from (7). 

At the next stage, we group together, over the whole lattice, the four factors 
with the same Ising spin (Tmn- These four factors come by the factorization of the 
four different bonds (weights) attached to a given mn site. Performing the averaging 
over amn ± 1 in each group, we find [7]: 

{ ^mn^mn-^rnnBmn } = 2 ^ (1 + ^Ln '^mn^^m-ln) (-'- "I" ^mn ^mn^mn-l) 
X(l + '^mn.O'mn) + ^rnnbran) = 1 + ^^mnP^mn '^m-lnPmn-l + '^mn^mn'^' (8) 

_LCf{i)„* _Lf(2)/,* \(ri -\-h ) 4- fW fC^) a * h* n h — 

— exp [ajnribmn + ^^vPmn^m-lriPmn-l + (^mn^m-ln + ^nirPmn-l) {^mn + ^mn)] ■ 

Taking also into account the diagonal Gaussian terms arising under factorization 
(7), we obtain the partition function (6) in the form [7]: 

L L L L 

Q ~ / n n dhmndbjnnd^mnd^mn ®^P { S S [^mn^mn bmn^mn'^ (n\ 
m=ln=l m=l n=l \'J J 



with the free-boundary condition for fermions: a^^ = 0, b^o = 0. Now Q is given as 
a Gaussian fermionic integral, equivalently, the 2D Ising model is reformulated as a 
free fermion field theory on a lattice. For space restriction, we do not comment here 
on the mirror-ordering procedure which just makes it possible to put the four relevant 
factors nearby as in (8), see for details [7,8]. It is important that the method works 
for the inhomogeneous lattices, this might be of interest, in particular, for studies 
of random systems [6,11]. The integral (9) can in fact be simplified integrating out 
the a^^, fields by means of the identity / dbda exp {ab + aL + L b) = exp (LL ), 
where L, L are some linear forms in Grassmann variables independent of a, b. 
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4. The 2D Ising models on complicated lattices. An important modifica- 
tion of the method has been introduced in [8] , where we start with the factorization 
of the cell weights presented by three-spin polynomials. The spin polynomial inter- 
pretation arises if we multiply few local weights forming elementary cell in Q. This 
enables us to obtain the fermionic representation for Q with only two fermionic vari- 
ables per site, which provides essential simplifications in the analysis [8]. For a set 
of 2D Ising models, including the standard rectangular, triangular and hexagonal 
lattices as the simplest cases, the partition function then reduces to the following 
expression [8]: 

Q — Sp { n (CKO + Oil (Tmn.O'm+ln + ^2 O'm+ln.O'm+ln+l + OI3 CTmn.O'm+l n+l) } , (m) 
(ct) mn \ / 

where ao, ai, q;2, oi^ are numerical parameters specific for each lattice (we now 
assume the homogeneous case). The three-spin polynomial in (10) is the effective 
Boltzmann weight of the elementary cell. The factorization method then yields the 
following Gaussian fermionic integral for the partition function [8]: 



<9 = / n dc;^Jc^^ exp { E [ ao CmnC„ - CmnC^-m - ^2 CmnC^n-l- .... 

mn mn ( J- J- ) 

~ 0^3 (^mnCjji-ln—l ~ "^l CmnCm—ln ~ ^mn^mn-1 \ S ' 



where c^„, is a set of the purely anticommuting Grassmann variables (two per 
cell). The derivation of (11) is simple [8]. 

The explicit evaluation of the integral (11) can be performed by passing to the 
momentum space (Fourier substitution for fermions). Taking then the limit of infi- 
nite lattice, we find the free energy per cell as follows [8]: 

lim^^^(^lng) = 

27r 27r J ^ 

= I / / f f In [ (a^ + a? + ai + ai) - 2 {a,a, - a^a^) cosp - (12) 
— 2 {aQa2 — CK^ag) cos q —2 {a^a^, — oc-^oc^ cos {p-\- q)] ■ 

The symmetries of this solution, and the closely related question on the location 
of the critical point, have an interesting interpretation within the spin-polynomial 
approach, as is discussed in [8], also see [12]. By a suitable specification of the pa- 
rameters aO)Cii)«2,a3 from (12) wc obtain, in particular, the solutions for the free 
energies for canonical rectangular, triangular, and hexagonal Ising lattices, as well 
as for some other lattices with more complicated local structures [8] . For application 
of these results to the analysis of the regularly diluted 2D Ising ferromagnets also see 
[9]. The fermionic integral for Q with the minimal number of the fermionic compo- 
nents (11) appears to be a suitable starting point to formulate the continuum-limit 
field theories for 2D Ising models near Tc, as is discussed in more detail in [12]. 
The resulting field theory is the massive two-component Majorana theory in the 2D 
euclidean space-time [12]. By doubhng of fermions, we can pass as well to the 2D 
Dirac field theory of charged fermions [12]. 

5. Conclusions. We have discussed some aspects of a simple non-combinatorial 
approach to the analytic solutions for the 2D Ising models on simple and more 
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complicated lattices. For any planar 2D Ising model the partition function can be 
expressed as a fermionic Gaussian integral. The integral over the anticommuting 
Grassmann fields is a powerful tool to analyze the two-dimensional Ising models. 
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